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Abstract 

The expression for the force on an ion in the presence of current can be derived from first principles without any assumption 
about its conservative character. However, energy functionals have been constructed that indicate that this force can be written 
as the derivative of a potential function. On the other hand, there exist compelling specific arguments that strongly suggest the 
contrary. We propose physical mechanisms that invalidate such arguments and demonstrate their existence with first-principles 
calculations. While our results do not constitute a formal resolution to the fundamental question of whether current-induced 
forces are conservative, they represent a substantial step forward in this direction. 
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Electrons flowing in a conductor can transfer momen- 
tum to the ions. From a static point of view, the steady- 
state electron density in a current-carrying conductor 
is microscopically different from that in the absence of 
current. In either picture, the passage of current gen- 
erates additional forces on ions. This may result in 
atomistic rearrangements and diffusion of atoms, a phe- 
nomenon known as electromigration. Q The advent of 
microelectronics has sparked off a considerable interest 
in electromigration since it has been found to be a ma- 
jor failure mechanism in metallic interconnects. In re- 
cent years, a new type of electronics is being considered 
which envisions the use of nanoscale components. Since 
nanoscale conductors can carry current densities much 
larger than their macroscopic counterparts, the problem 
of electromigration is again of key importance. Recent 
work has studied the effect of current-induced forces in 
atomic S S H E S S H and molecular wires [13, El- 
Combined with early investigations on the subject, [13 
this work has shed new light on the microscopic mecha- 
nism of current-induced mechanical effects. 

While a consensus on the physical origin of current- 
induced forces seems to have been reached, a fundamen- 
tal question remains unanswered: are current-induced 
forces conservative? The difficulty in answering this ques- 
tion is intimately related to the difficulty in defining a 
force in a non-equilibrium problem, [l^ [Q 

In this paper we address this question in the follow- 
ing way: we consider representative arguments against 
the conservative nature of current-induced forces and 
suggest physical mechanisms that invalidate these argu- 
ments. We then corroborate the existence of such mech- 
anisms with first-principles calculations. While this by 
no means proves that these forces are conservative, we 
believe it is a substantial step forward in this direction 
and hope that our suggestions will generate much-needed 
discussions. 

Let us start from the expression for forces in the pres- 
ence of current. This expression can be derived in dif- 
ferent ways. One approach is to start from a fully dy- 



namical quantum description of the many-body problem 
comprising electrons and ions, and then take the limit as 
the ionic motion becomes classical. 13j In this derivation, 
the force expression comes about as the classical limit 
of Ehrenfest's theorem applied to the rate of change of 
the momentum of an ion. 13] Formally, this limit cor- 
responds to high ionic masses or high kinetic energies 
[iq . A different approach is to start from an appropriate 
mixed quantum-classical Lagrangian. The force expres- 
sion then results from the Euler-Lagrange equation of 
motion for the classical ions . Both approaches result 
in the same expression for the force on a given ion, with 
position R, due to the self-consistent electronic density 
p(r) under current flow 16], 



which can be equivalently written as [l^ 
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where v is the electron-ion interaction. The first term 
on the RHS of Eq. Q is similar to the usual Hellmann- 
Feynman contribution to the force due to localized elec- 
tronic states \tpi). The second term is the contribution 
to the force due to continuum states. 13] The wavefunc- 
tions |?/^a) are eigendifferentials in the continuum a. 
Due to the algebraic equivalence of the above two equa- 
tions, in the following, when discussing current-induced 
forces, we will refer to either one without distinction. 

Let us now turn to the question of whether or not 
these forces are conservative. It is essential to emphasize 
that the force obtained either from the Ehrenfest's the- 
orem 13] or a Lagrangian formulation 16], is not neces- 
sarily conservative because the above first-principles dy- 
namical derivations do not contain any assumption about 
the existence, or otherwise, of a free energy functional, 
of which this force might be the derivative. However, as 
a separate exercise, it was suggested both in Ref. 13 and 
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FIG. 1: Schematic of a current- carrying wire between two 
bulk electrodes. An atom is moved from a point A to a point 
B along two different paths indicated as solid lines, d is the 
largest distance of the atom from the wire along the vacuum 
path. Electrons flow from right to left. 



in Ref. ^ that forces in the presence of current can also 
be written as derivatives of appropriate energy function- 
als. For instance, in Ref. fl3 the stationarity of an action 
was used to map the force problem into a variational 
one in which the force can be written as the derivative 
of an energy function. In Ref. 14 an entropy function 
was used, whose maximization (constrained to the en- 
ergy of the current-carrying system and the numbers of 
electrons in it that originate from the right and left reser- 
voirs) again yields the force as the derivative of a poten- 
tial function. These independent derivations suggest that 
current-induced forces are conservative. 

On the other hand, there exist compelling specific 
arguments that strongly suggest that current-induced 
forces cannot be conservative. These arguments consist 
in gedanken experiments, aimed at showing that the work 
done by current-induced forces along specific closed paths 
is not zero and that, therefore, these forces are not conser- 
vative. Below we refer to one notable such construction 
due to Sorbello. ^ This example is of charming simplicity 
and is particularly convincing. The argument goes like 
this. Picture a single defect atom in a current-carrying 
wire between bulk electrodes under an applied bias 
V . First, let the atom move, quasi statically, from a point 
A to a point B inside the wire (see schematic in Fig.^. 
The current-induced force on the defect in the wire does 
some non-zero work on the way. Now imagine taking the 
atom from Aio B along a path that goes outside the wire 
(see Fig. Since outside the wire there apparently is 
no current, no current-induced force should be present. 
Hence, the work done by current-induced forces on the 
atom along the path from Aio B inside the wire cannot 
be the same as along the path outside the wire. This 
would be sufficient to prove that current-induced forces 
are not conservative. 

However, there is something missing in this argument: 
the quantum-mechanical nature of current outside the 
wire. To see this, picture the atom at a point in space 



outside the wire. Let be the global current-carrying 
electron states for the system. If there always is some 
(however small) quantum-mechanical coupling between 
the atom and the wire, there will always be some non- 
zero tunneling probability for electrons to get from the 
wire onto the atom, and back. The atom then acts as 
a resonant center and there will always exist an energy 
window where the states contain resonant paths that 
emerge from the wire, pass through the atom, and go 
back into the wire. Due to the scattering of electrons at 
the atom, these resonant paths may be expected to cre- 
ate a charge dipole, centered at the atom, much like the 
local resistivity dipole in electromigration. 1] Increas- 
ing the distance d between the atom and the wire will 
make the resonant energy window narrower, while the 
resonance will remain pinned against the energy window 
eV, due to self-consistency. Now, the total current car- 
ried by the resonance is given by the energy integral, 
over the energy window eV ^ of a transmission function, 
T{E). As the resonance gets narrower in energy, so does 
T{E). T{E) is always bounded and hence, as d increases, 
the total current carried by the resonance will decrease. 
On the other hand, the charge dipole across the resonant 
center originates from the density of states, DOS(£^), due 
to the populated current-carrying electron states in the 
energy window eV. By contrast with T{E)^ DOS(^) is 
unbounded: if d increases, the resonant peak in DOS(£^) 
gets narrower but taller in such a way that its energy 
integral remains constant. Since current-induced forces 
are related to the integral of the density of states [see 
e.g. Eq. we expect the force on the atom due to the 
resonant dipole to be constant with d. In other words, 
there is always a current-induced force on the atom, even 
outside the wire, and this force does not vanish with in- 
creasing d. 

A key assumption of our reasoning is that there always 
is some coupling between the atom in the vacuum and 
the metal. If this coupling is truncated artificially, then 
the current and the "resonant-dipole force" will both be 
discontinuously extinguished. However, such a trunca- 
tion would at the same time create a bound state (at the 
atom) in the continuum, and would thus modify the na- 
ture of the electronic spectrum and of the corresponding 
Born-Oppenheimer (BO) manifold. Such a discontinu- 
ous modification corresponds to a breakdown of the BO 
approximation, in which case the notion of conservative 
forces is lost by definition. 

Our analysis above relies only on generic properties of 
resonant transmission. We therefore now seek to demon- 
strate the existence of a resonant dipole in the simplest 
possible resonant configuration: that of a single atom in 
the vacuum between two bulk electrodes. Our goal is 
to show that this resonant dipole gives rise to a finite 
current-induced force that does not vanish with increas- 
ing electrode separation, even when the current dies out. 
To this end, we have performed first-principles calcula- 
tions of the current-induced force on a single Si atom 
between two jellium electrodes. The interior electron 
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FIG. 2: Density of states due to the populated global current- 
carrying states originating from the right electrode for the 
system with the Si atom, minus the same quantity for the 
electrode-electrode system without the atom, at 0.1 V and for 
different separations of the jellium edges. Zero of energy cor- 
responds to the right electrochemical potential. The vertical 
dotted line indicates the position of the left electrochemical 
potential. 



density of the electrodes is taken equal to the value for 
metallic gold (r^ ~ 3). 20] The electrodes are kept at a 
certain bias with the left electrode positively biased so 
that electrons flow from right to left. In Fig. [2 we plot 
the density of states as a function of energy for a bias of 
0.1 V and different distances between the jellium edges of 
the two electrodes. It is clear from this figure that a res- 
onance develops in the energy window between the left 
and right chemical potentials. The resonance becomes 
narrower and taller with increasing distance. It is also 
evident from Fig.j^lthat this resonance is pinned towards 
the right chemical potential. Such pinning can be intu- 
itively understood as follows. By increasing the distance 
between the electrodes the resonance becomes narrower 
and narrower. In order to preserve charge neutrality, the 
amount of charge accommodated by the resonance needs 
to be the same at all distances. The only way to achieve 
this is for the "center" of the resonance to move closer to 
the right electrochemical potential. 

Now that we have shown the existence of the reso- 
nance, let us look at the spatial distribution of the charge 
dipole it creates across the atom. In order to do this we 
follow Lang's approach 3] and calculate first the differ- 
ence between the electronic distribution at a given bias 
p(r, V) and the corresponding distribution at zero bias 
p(r,0). Call this quantity p(r, V). It contains three con- 
tributions: (i) the polarization of the biased electrode- 
electrode system without the atom, (ii) the polarization 
of the electron cloud at the atom due to the field gener- 
ated by (i), and (iii) the current-induced resonant dipole. 
The first two terms generate a force of equal magnitude 
but opposite sign so they do not contribute to the total 
force. (An isolated atom in a uniform electric field does 
not have a net force on its nucleus). To eliminate (i) 
and (ii), we subtract (i) from p(r, F) and take the dif- 
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FIG. 3: Current-induced resonant vacuum dipole as defined 
in the text for a jellium- jellium edge distance of 10 a.u. Dark 
area on the right of the nucleus (represented by a dot) corre- 
sponds to negative charge. 



ference between the result of that subtraction at a given 
electrode-electrode distance and the same quantity at a 
larger separation scaled to the field of the previous dis- 
tance. 3] This gives a negative fraction of the current- 
induced dipole, so a change of sign is necessary. Call this 
quantity p/(r, V). Finally, to isolate the antisymmet- 
ric (dipole) component of this distribution we evaluate 
[p/(r, V) - p/(r, -V)]/2. The result is plotted in Fig. H 
for a separation of 10 a.u. (with respect to 16 a.u.) and 
clearly shows the current-induced resonant dipole across 
the atom. The dipole has positive charge on the left of 
the atom and negative charge on its right so that the 
current-induced force pushes the atom to the right, i.e. 
against the electron flow. 

Finally, we have to show that the current decreases 
with electrode separation while the force does not. This 
is plotted in Fig.^for different electrode separations. It is 
clear from Fig. ^ that the current decreases substantially 
with electrode separation while the force saturates at a 
value of about 2 mRy/a.u. Notice again that this force is 
positive, i.e. it pushes the atom against the electron flow, 
in agreement with the sign of the resonant dipole. These 
numerical results confirm our conjecture: as long as the 
atom is coupled to the electrodes, a resonant dipole forms 
across it, such that a non-zero current-induced force is 
exerted on the atom that persists with separation even 
as the current dies out. 

We conclude by stressing again that our goal here is 
to show that existing arguments against the conserva- 
tiveness of current-induced forces ignore specific physical 
mechanisms. These mechanisms predict the existence of 
current-induced forces in unexpected places, such as an 
atom in the vacuum next to a current-carrying metal. 
We have demonstrated the existence of such forces using 
first-principles calculations. These results do not prove 
that current-induced forces are conservative but we be- 
lieve that they give extra support in favor of such a con- 
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FIG. 4: Current (right scale) and current-induced force (left 
scale) for a single Si atom between two bulk electrodes as a 
function of distance between them. 



elusion. More theoretieal work is elearly neeessary to 
fully resolve this fundamental problem. 
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